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Abstract 

Nucleation in the two-dimensional (;-state Potts model has been studied by 
means of Monte-Carlo simulations using the heat-bath dynamics. The initial metastable 
state has been prepared by magnetic quench of the ordered low-temperature phase. 
The magnetic field dependence of the nucleation time has been measured as the 
function of the magnetic field for different q and lattice sizes at T = 0.5 T^. A size- 
dependent crossover from the coalescence to nucleation region is observed at all q. 
The magnetic field dependence of the nucleation time is roughly described by the 
classical nucleation theory. Our data show increase of the anisotropy in the shape 
of the critical droplets with increase of q. 
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1 Introduction 

Nucleation of 'bubbles' or 'droplets' of the stable phase is a very common phenomenon 
which initiates relaxation of metastable states near the first order phase transition in 
many systems of condensed matter physics quantum fields 0], and cosmology P, ^. 
The main quantity of interest for the theory and applications is the nucleation time, which 
characterizes the time required by the system to leave the metastable state. In lack of 
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a consistent theory of the metastabihty, the nucleation problem has been subjected to 
extensive study in computer simulations. Numerical data are usually interpreted in terms 
of the phenomenological classical nucleation theory, or in terms of its field-theoretical 
version proposed by Langer |^. A great deal of this work has been carried out on the 
Ising model supplied with some appropriate dynamics. 

It seems reasonable, however, to extend the numerical study of the nucleation to 
more complex and rich models, taking into account a great variety of physical realizations 
of this phenomenon. In this paper, we have studied nucleation in the two-dimensional 
g-state Potts model by use of the Monte-Carlo simulation. Due to a rich equilibrium 
phase structure, this model provides a lot of possibilities to form the initial metastable 
phase. For g > 4, the zero-field Potts model has the first order phase transition point 
driven by the temperature 0. Near this point, the metastable state can be prepared by a 
quenched cooling from the high to low temperature phase. Relaxation of such a state in 
different regimes has been studied numerically by Safran et. a/ ^, and by Arkin et. al 
T0| , |lT[. Meunier and Morel analysed metastabihty in the thermally quenched Potts 



model within the Fisher's droplet picture approach Arkin et. al confirmed UTO], [IT 



that in a case of slight quench (small |Tc — T|), the equilibrium proceeds through the 
usual nucleation mechanism. However, lack of the macroscopic order parameter (like the 
average magnetization in the Ising system below Tc) in the both initial and final states 
hinders precise measurement of the nucleation time in this regime. 

We simulate relaxation of the metastable states, which are induced in the Potts fer- 
romagnet near the first order phase transition driven by the external magnetic field H. 
The heat-bath dynamics is used to thermalize spins. Nucleation time is measured as the 
function of the applied field at T = 0.5 Tc for g = 2, 3, 5, 15. 

One can say, that for small g , g < 4, the line H = 0, < T < in the Potts model 
represents the liquid-vapour phase transition, and it simulates the liquid-solid transition 
for large g, g > 4. It should be noted, that the hanger's theory of nucleation based on 



the coarse-grained free-energy functional can not be applied in the latter case So, 
our aim is to identify in the Potts model the nucleation regimes well-known in the Ising 
model, and to study how they are modified by increase of g. 

2 The Potts Model 

We study the g-state ferromagnetic Potts model on the square lattice having L rows and 
L columns. The discrete spin variable a = 1, 2, g is associated with each lattice site j. 
Helical boundary conditions are used [Q. The model Hamiltonian is defined as 

<ij> j 
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Here the first summation is over the nearest neighbour pairs, 2J > is the couphng 
constant, H is the external magnetic field applied along the 2-direction, and 



1 ifa = a', 
lia^a'. 



We use also dimensionless parameters K = J (3, and h = Hf3, where /3 = l/fc^T is the 
inverse temperature. The macroscopic state of the lattice can be described by q real 
positive order parameters M{a) 

M{a)=L-^Y.^{a,,a), (2) 
j 

cr=l 

The model undergoes ferromagnetic phase transition at the critical temperature 

2J 



(3) 



The ferromagnetic low-temperature phase at zero field is g-times degenerated. Different 
ferromagnetic phases are enumerated by the discrete parameter a. This implies that in 
the (j-phase, the orientation a dominates in the lattice. For a review of many others 



known properties of the Potts model see [15 



3 Classical Nucleation Theory 

The classical nucleation theory (CNT) considers thermally induced generation of the 
stable phase droplets in the metastable surrounding. A droplet of radius R can appear 
with probability p{R) proportional to the Boltzmann factor 

p{R) ~ exp[-/5^(i?)], 

where T{R) is the droplet free energy^: 

T{R) = 2rrR s - irR^M H. (4) 

Here s is the surface tension on the droplet boundary, M is the zero-field spontaneous 
magnetization. The critical droplet radious Rc maximizes the free energy @. It is postu- 
lated in CNT, that the nucleation time t{h, L) is inversely proportional to the probability 
of generation of the critical droplet: 

t{h,L) ~ L-^ exp[(3T{R,)] = L'^ exp[7r{l3sf/{M h)]. (5) 

^In the Ising model, the volume free energy —2ttB?M H of the droplet has an additional factor 2 due 
to a different definition of the magnetic field H . 
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The classical nucleation theory considers only spherical droplets. This is reasonable, if 
the surface tension of the domain wall does not depend on its is orientation. If the surface 
tension is anisotropic, the critical droplet are supposed to have the equilibrium shape 
given by the Wulff's construction 0]. In this case equation (|) should be modified to 

~ L-2exp[A(T)//i], (6) 

where the amplitude A{T) is determined by the equilibrium droplet shape. In the Ising 
model, this shape is given by the equation ITT 



coshxi + cosha;2 = sinh(2_ft') H , (7) 

smh(2ii ) 

where xi, X2 denote Descartes coordinates on the droplet boundary. It varies from circle 
at T = Tc — to square at T = 0. 

The amplitude A(T) can be written as A(T) = S(T)/M, where S(T) denotes the area 
bounded by the curve (^, and 

M = {1- [smh{2K)]'^y/\ 
Numerically we have for the q = 2 state Potts (Ising) model at T = 0.5 T^: 

5(0.5 = 6.537, M(0.5TJ = 0.998, 

A(0.5TJ = 6.55. (8) 

It should be noted, that equation (D gives only the leading /i — > asymptotics of the 
nucleation time. The most important correction is expected to have the form of the 
prefactor with a power-law /i-dependence 

t{h,L)=L-'^exp[A{T)/h]. (9) 
The exponent b here arises in the field theoretical nucleation theory from the contribution 



of the capillary waves on the critical droplet surface |T8[, 6 = 1 in the two-dimensional 
Ising model. The exponent c gives the h variation of the so-called kinetic prefactor, 
which depends on the chosen dynamics. For dynamics which can be described by the 
Fokker-Plank equation, the value c = 2 is expected [0, |20[| . 

The exact equilibrium (Wulff's) shape of the droplet is not known for the Potts model 
with g > 3. So, we had to use more simple and less accurate estimates for the amplitude 
A{T). We have used square and circle approximations for the critical droplet shape. The 
surface tension in the both cases was replaced by its zero-temperature value s = 2 J. This 
yields for the amplitude A{T) in the nucleation time exponent: 



16 K"^, for square droplet, 
Att . for circle droplet 



Table 1: Amplitude A(T) at T = 0.5 Tc for square and circle droplets. 



q 


A(0.5T,) 




square 


circle 


3 


16.2 


12.7 


5 


22.1 


17.3 


15 


40.1 


31.5 



Numerical values for the amplitude A{T) calculated from (|T0|) for different g at T = 
0.5 Tc are given in Table |l[ 

The size-dependent behaviour (^) of the nucleation time is observed in not very large 
lattices, in the so-called nucleation regime. In this regime, only one droplet of the stable 
phase appears in the system, and then expands to the whole lattice. 

At large L, the coalescence regime is realized. In this regime the relaxation time does 
not depend on L. Many critical droplets appear in the system, expand and coalesce. 

It is expected from the Kolmogorow- Johnson- Mehl-Avrami theory |jl6|, that 

t{H)^exp[A{T)/{3h)] (11) 
in the coalescence regime in two-dimensional systems . 

4 Simulation Scheme 

In the Potts model, the applied field and the initial metastable configuration of spins in 
the lattice can be chosen in many different ways. We have studied only two of them. 

• First regime. A positive magnetic field if > is applied along the 2-direction, 
whereas all the spins are initially oriented along the 1-directions: cxj = 1. Such a 
field removes completely degeneracy of the low temperature phase, leaving only one 
equilibrium ferromagnetic state: cr = 2. 

• Second regime. A negative magnetic field if < is applied along 2-direction, and 
all spins are initially oriented in the same direction aj =2. The ferromagnetic phase 
still remains {q — l)-times degenerated (if g > 3) in such a field. In this regime, the 
system moves to some bifurcation point, and then relaxes to one of the equivalent 
(g — 1) equilibrium states. 

The heat-bath technique (Glauber kinetics) has been used to orient the spins. To produce 
a sequence of pseudorandom numbers, we have used the standard generator of pseudo- 
random numbers with multiplication by 16807 to give 64-bit integers. The candidate 
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spins to flip have been chosen sequentially: site by site in the row, and row by row in 
the lattice. We define the nucleation time t as the number of sweeps through the lattice 
after which the number of spins in the initial (metastable) phase reduces to 0.5L^. The 
/i-dependence of the nucleation time t has been measured at L = 31, and L = 1001. In 
the nucleation regime, the nucleation time fiuctuates very much. So, we had to average 
the nucleation time in this regime over 100 or even 1000 runs. It is important to note, 
that for q > 3 the relation between the nucleation time t and the total relaxation time to 
is not straightforward and is different in different regimes. 

A typical time dependence of the order parameters M{a) in the first nucleation regime 
for g = 3 is shown in Fig. 1. Since the fall of the magnetization M(l) is very rapid 
compared with the time which the system has spent in the metastable state, the nucleation 
time t and the total relaxation time to are essentially the same in this regime. This is 
true also in the second nucleation regime. In the latter case, however, the system can fall 
with equal probabilities to one of the equilibrium states, either a = 1, or a = 3. 

Fig. 2 shows a typical time evolution of the magnetization in the first coalescence 
regime. Up to the saturation at 1, the magnetization M(2) increases approximately 
linearly with time due to the steady expansion of many stable phase droplets. Thus, the 
system reaches equilibrium in about two nucleation times: to — 2 1. 

The situation is quite different in the second coalescence regime, see Fig. 3. In this 
case many critical droplets of the both equilibrium phases a = and a = 1 are initially 
generated in the lattice. They expand and expel the metastable a = 2 phase from the 
system in about two nucleation times. The resulting state is still far from the equilibrium, 
since the lattice is divided into domains of two stable phases a = 1 and a = 3. During the 
final (very slow at low temperatures) stage of the relaxation, these two phases compete, 
until one of them will expel another from the system. So, the nucleation time t and the 
total relaxation time to may be of different order of magnitude in the second coalescence 
regime. 



5 Results 

Our simulations were restricted to the temperature T = 0.5 and the lattice sizes L = 31 
and L = 1001. We have started from the q = 2 Potts model, which is equivalent to the 
Ising model. This provides a good possibility to verify efficiency of the applied simulation 
scheme, since nucleation in the Ising model has been thoroughly studied JI], |23|, |2^ . 



Figure 4 shows the nucleation time in the logscale for the q = 2 Potts (Ising) model. 
Pluses correspond to the nucleation time averaged over 1000 runs in the L = 31 lattice. 
Stars relate to the nucleation time in the L = 1001 lattice. Crosses indicate the dispersion 
D[t] of the nucleation time distribution in the L = 31 lattice: 



where (...) denotes averaging over different runs. 

Three relaxation regimes are clearly identified: the strong field regime (SFR), the 
coalescence regime (CR), and the nucleation regime (NR). Increase of the lattice size does 
not change the nucleation time in the strong field and coalescence regimes. The nucleation 
regime was not observed in the L = 1001 lattice in the studied region of fields h. In the 
nucleation regime, the nucleation time fluctuates very much, and the distribution p{t) of 
these fluctuations is described by the simple exponential law: 

p(t) = -exp(-t/r). (12) 
r 

For this distribution, the averaged time equals to the dispersion: {t) = D[t] = r. This is 
indeed observed in the L = 31 lattice at 1/h > 1.15 indicating well established nucleation 
regime in this fleld region. The theoretical prediction for the averaged nucleation time 
in the nucleation regime with the Wulff's shape droplet is given by equation (|^), where 
6 + c = 3, and A = 6.55. The upper straight line in Fig. 4 shows the nucleation 
theory exponent, and the second (curved) line plots the the same exponent corrected by 
the prefactor Ch~^, i.e. equation @. Though the latter curve contains only one free 
parameter - the constant factor C - it flts well the numerical simulation points for the 
L = 31 lattice in the nucleation regime, much better then the flrst (straight) line. 

The curved line gives also a good flt for the dispersion points in the small L = 31 
lattice not only in the nucleation region, but also at crossover to the coalescence regime 
down to ^ 0.8. The possible explanation of this result may be the following. Measured 
nucleation time t consists of two terms t = tn+te, which have the same order of magnitude 
in the crossover region. The flrst term relates to the time of nucleation of the flrst critical 
droplet in the system, and the second term denotes the time of the expansion of this 
droplet to 1/2 of the whole lattice. The flrst term should strongly fluctuate according to 
the exponential law (12). If the second term te fluctuates much weaker, one should have 



for the total time dispersion: 

D[t\ =D[tn + te] ^ Dm = (t„) = r, 

where the averaged time r of the flrst droplet nucleation is still given by equation (^. 
Thus, in the crossover region, dependence (|]) should flt better the time dispersion points 
D[t], than the averaged time (t), which is indeed observed. 

The flrst numerical veriflcation of the magnetic fleld dependence of the prefactor in 
(|^) for the two-dimensional Ising model has been reported by Rikvold et. al ||2^. In this 



work, the value 6 + c = 3 was conflrmed at T = 0.8 for the Metropolis dynamics with 
randomly chosen sites. However, a smaller value 6 + c ~ 2 was observed in , when 
the sites in the lattice were chosen sequentially. Though we have also used the sequential 
chose of sites in the lattice, our data for the nucleation time in the nucleation region are 
better fltted by (0) with the prefactor exponent 6 + c = 3, rather then 6 + c = 2, see Fig. 



5. This is more clear for the dispersion D\i\ (crosses in Fig. 5), since the latter is still 



The lower straight line is the best linear fit for the nucleation time in the larger 
L = 1001 lattice in the coalescence regime (CR). The ratio of slopes R = 1.9 of two 
straight lines in Fig. 4 corresponding to the nucleation and coalescence regimes is much 
smaller than 3, expected from the theory in the small h limit. However, we are still far 
from the well-established small h limit, as it is evident from the considerable difference 
between first and second lines in Fig. 4. If one takes this into account, this disagreement of 
the observed R = 1.9 value with 3 expected in the h ^ limit becomes not so surprising. 

Figures 6-10 show the observed dependence of the nucleation time on the inverse field 
in the g = 3, 5, 15 Potts models, L = 31, 1001 in the first h > and second h < 
regimes. Pluses correspond to the nucleation time averaged over 100 runs in the smaller 
lattice. Crosses indicate nucleation time in the L = 1001 lattice. Only the second regime 
h < has been studied in the g = 15 Potts model, since the nucleation time in the first 
regime h > becomes too long in this case, see equation (|13|) below. 

Qualitatively, the magnetic field dependence of the nucleation time in all cases is 
similar to that in the Ising model. The same three field regions (SRF, NR, CR) have been 
observed. As in the Ising case, the nucleation regime was realized only in the smaller 
lattice, L = 31. In the latter regime, the relaxation time strongly fluctuates. 

In all cases, log(t) behaves linearly with l//i in the nucleation regime (NR), in agree- 
ment with the classical nucleation theory. In the g = 3 and g = 5 systems the observed 
slopes (log(t) vs. h^^ ) in NR agree with the circle shape of the critical droplets. However, 
in the g = 15 model our data indicate the square critical droplet shape, see Fig. 10 , and 
Table |l|. 

The ratio of slopes in nucleation and coalescence regimes varies from R = 2.6 in the 
q = 3, h > case, down to R = 2.1 in the q = 5, h < case. As in the g = 2 model, these 
values are lower than 3 predicted by the theory. 

In the first nucleation regime {h > 0), the nucleating critical droplet should be of 
the single equilibrium phase a = 2, whereas in the second nucleation regime {h < 0), a 
critical droplet of either of (g — 1) equilibrium phases a, a ^ 2 provides relaxation. So, 
it is naturally to expect, that nucleation time (ti) in the first nucleation regime will be 
longer, that the nucleation time (^2) in the second nucleation regime, and 



The observed ratios (ti)/ (^2) determined from the circle droplet fits in Fig. 6-9 



described by (|^) at the crossover to the coalescence regime. 




exp(-2.95)/ exp(-3.64) = 2.0 for g = 3, 
exp(-1.29)/ exp(-2.8) = 4.5 for g = 5 



(14) 



are in a good agreement with (|T3|). 
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6 Summary 



We have simulated nucleation in the g-state Potts model on the square lattice for q = 
2, 3, 5, 15. The magnetic field dependence of the nucleation time have been measured in 
two regimes {h > and h < 0) at T = 0.5 T,, for two lattice sizes L = 31 and L = 1001. 
The strong field regime, coalescence regime, and nucleation regime are clearly identified. 

The observed magnetic field dependencies of the nucleation time is in agreement with 
the classical nucleation theory, if one takes into account the anisotropy of the critical 
droplet shape. Our numerical data indicate, that the shape of the critical droplet changes 
from (roughly) circle towards square with increasing q at the fixed temperature T = 



0.5 Tc. A similar result was observed by Bikker et. al in numerical calculations of 
the equilibrium crystal shape in the three-dimensional Potts model. So, increase of q has 
a similar effect on the equilibrium droplet shape, as decrease of temperature. This is 
natural, since the both reduce the correlation length, which should increase anisotropy in 
the system. 

The ratio of (\og{t) vs. h^^) slopes in nucleation and coalescence regimes does not 
approach 3, predicted by the theory in the h —>■ limit. We associate this discrepancy 
with the 'finite-/;, corrections', which are still considerable in the studied region of the 
magnetic field. 

Though g = 2, 3 and g = 5, 15 state Potts models have different type phase transitions 
at Tc (continuous for g = 2, 3 and first order for g = 5, 15), metastable phases in these 
models relax to the equilibrium in a similar way. 
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Figure 1: Decay of metastable state in first nucleation regime 
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Figure 2: Decay of metastable state in first coalescence regime 
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Figure 3: Decay of metastable state in second coalescence regime 
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Figure 4: Nucleation time and dispersion plotted against 1/h va. the 
q — 2 Potts (Ising) model for L = 31 
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Figure 5: Averaged nucleation time {t) and dispersion D[t\ divided by 
exp[6.55//i], versus 1/h (loglog-scale) in the q — 2 Potts (Ising) model. 
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Figure 6: Nucleation time versus 1/h in the q = 3 Potts model, h > 0. 
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Figure 7: Nucleation time plotted against l/\h\ in the q — 3 Potts 
model, h < 0. 
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Figure 8: Nucleation time plotted against 1/hm the q — 5 Potts model, 
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Figure 9: Nucleation time plotted against l/\h\ in the q 
model, /i < 0. 
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Figure 10: Nucleation time plotted against l/\h\ in the q — 15 Potts 
model, h < 0. 



15 



